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Abstract. We present the first systematic spectropolarimetric study of Luminous Blue Variables (LB Vs) in the Galaxy and the 
Magellanic Clouds, in order to investigate the geometries of their winds. We find that at least half of our sample show changes 
in polarization across the strong Ho- emission line, indicating that the light from the stars is intrinsically polarized and therefore 
that asphericity already exists at the base of the wind. Multi-epoch spectropolarimetry on four targets reveals variability in their 
intrinsic polarization. Three of these, AG Car, HR Car and P Cyg, show a position angle (PA) of polarization which appears 
random with time. Such behaviour can be explained by the presence of strong wind-inhomogeneities, or 'clumps' within the 
wind. Only one star, R 127, shows variability at a constant PA, and hence evidence for axi-symmetry as well as dumpiness. 
However, if viewed at low inclination, and at limited temporal sampling, such a wind would produce a seemingly random 
polarization of the type observed in the other three stars. Time-resolved spectropolarimetric monitoring of LBVs is therefore 
required to determine if LBV winds are axi-symmetric in general. 

The high fraction of LBVs (> 50%) showing intrinsic polarization is to be compared with the lower ~ 20-25 % for similar 
studies of their evolutionary neighbours, O supergiants and Wolf-Rayet stars. We anticipate that this higher incidence is due to 
the lower effective gravities of the LBVs, coupled with their variable temperatures within the bi-stability jump regime. This is 
also consistent with the higher incidence of wind asphericity that we find in LBVs with strong Ha emission and recent (last ~ 
10 years) strong variability. 

Key words. Techniques: polarimetric - Stars: mass-loss - Stars: winds, outflows - Stars: early-type - Stars: activity - Stars: 
evolution 



1. Introduction 

Luminous Blue Variables (LBVs), or S Dor variables, are a 
class of evolved massive stars located close to the empiri- 
cal upper-luminosity boundary on the H-R diagram, known as 
the Humphreys-Davidson (HD) limit (Humphreys & Davidson 
1994; van Genderen 2001, hereafter vGOl). They exhibit pho- 
tometric variability on three different scales: microvariations 
of a few tenths of a magnitude on timescales of weeks to 
months; variations of ~ 1 mag on timescales of years; and 
massive eruptions of > 2 mag, observed in only a handful of 
objects, rj Car and P Cyg being the most famous (see Lamers 
et al. 1998; Humphreys et al. 1999). At visual minimum their 
spectra resemble those of early-type supergiants with emission 
lines of H i, He i and Fe u, with temperatures ranging from 
15-30,000K. At visual maximum, their temperatures are ob- 
served to fall to ~ 9,000K, as their radius increases and atmo- 
sphere cools (de Koter et al. 1996). Their mass-loss rates, of ~ 
10~ 5 - 10~ 4 M yr~', vary with effective temperature, which can 
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be understood by radiation pressure on spectral lines (Vink & 
de Koter 2002). The LBV phase itself is expected to last ~ 10 5 
yrs, and is thought to be a transitional phase experienced by 
stars with an initial mass > 25M prior to the Wolf-Rayet stage 
(e.g. Maeder 1997). 

With only a handful of exceptions, LBVs are surrounded by 
expanding aspherical nebulae (Weis 2003; Nota et al. 1995). 
These nebulae can be bipolar, elliptical or irregular; be dis- 
crete shells or clumpy, and are thought to result from the stars' 
episodic mass-loss phases. Nota et al. (1995) proposed that 
these nebulae may be formed by a spherically symmetric wind 
interacting with a pre-existing density contrast left over from 
a prior mass-loss episode, which was subsequently backed-up 
with models (Frank et al. 1995; Dwarkadas & Balick 1998; 
Langer et al. 1999). It has however also been shown that such 
nebulae can be formed if the wind itself is axi-symmetric. 
A bipolar wind can be created if the rotation rate is close to 
break-up and the star's flux is latitude dependent due to oblate- 
ness and gravity-darkening (e.g. Dwarkadas & Owocki 2002). 
An equatorial wind may exist if the star is close to the bi- 
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Table 1. Log of the LBV observations. The signal-to-noise (S/N) was calculated from the noise in flat regions of the continuum. The spectral 
resolution was calculated using the FWHM of lines in the arc calibration spectra taken at the time and position of the observed star. All 
observations, including the archive data, were done at the AAT except P Cyg which was done at the WHT All data was taken using a 1200R 
grating, except HD 160529 on 18/9/02 which was done using a 600V grating. 



stability jump temperature where the recombination of Fe iv 
to Fe m below the sonic point of the wind leads to an increased 
opacity and a marked change in wind properties (Vink et al. 
1999). Indeed, such a scenario has been suggested by Lamers 
& Pauldrach (1991) and Pelupessy et al. (2000) to explain the 
two-component wind model of the possibly related B[e] su- 
pergiants. However, direct observational investigations into the 
present-day wind geometries of evolved massive stars in gen- 
eral, and LBVs in particular, are thus far thin on the ground. 

Evidence for aspherical stellar winds can be found through 
spectropolarimetry. Ionised circumstellar gas in the wind which 
has a flattened geometry on the sky, such as in a bi-polar or 
equatorially-enhanced flow, (electron) scatters continuum pho- 
tons originating from the star. In the optically thin case, the 
scattered continuum light is then linearly polarized perpendic- 
ular to the plane of the flow. In the case of an optically thick 
wind, multiple scattering effects can mean the light is polar- 
ized parallel to the plane of the flow (Angel 1969; Wood et al. 
1996). Emission-line radiation emanating from the ionized gas, 
which is formed over a much larger volume at a larger radius, 
undergoes less scattering and remains essentially unpolarized. 
Therefore a drop in polarization across an emission line is in- 
dicative of aspherical structures within the line-forming region. 
As the bulk of the polarization occurs within a couple of stel- 
lar radii where the circumstellar gas is most dense (Cassinelli 



et al. 1987), this signature tells us about the geometry at the 
very base of the wind. 

Spectropolarimetric evidence for aspherical stellar winds 
exists for the Galactic LBVs AG Car (Leitherer et al. 1994; 
Schulte-Ladbeck et al. 1994), HR Car (Clampin et al. 1995), 
and P Cyg (Taylor et al. 1991; Nordsieck et al. 2001), as well 
as the Magellanic Cloud LBV R 127 (Schulte-Ladbeck et al. 
1993). These studies deal only with individual objects, how- 
ever, and until now no systematic study exists to determine 
whether LBVs in general undergo aspherical mass loss. Such 
studies do exist for O supergiants and Wolf-Rayet stars - of 
which LBVs are possibly the evolutionary mid-point - but ev- 
idence for significant intrinsic polarization was found in only 
~ 25% and ~ 20% of these groups respectively (Harries et al. 
2002, 1998). As LBVs are closer to their modified Eddington 
limit, have lower effective gravities, and are in the bi-stability 
jump regime, we may expect to find different mass-loss be- 
haviour. In the first systematic spectropolarimetric study of 
LBVs, we will show that they do indeed show a higher propor- 
tion of intrinsic polarization, and hence that their wind geome- 
tries in general differ from those of their evolutionary neigh- 
bours. 

We begin in Section 2 by describing the observational tech- 
nique and data-reduction steps. We describe in Section 3 an 
empirical method for estimating the stars' temperatures (and 
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Table 2. Log of the observations of the stars in the spectral atlas. Columns show the same as Table 1, except column 4 which shows the MK 
spectral type of each of the stars. All observations were done at the AAT, and all used a 1200R grating with the exception of HD 183143 which 
used a 600V grating. All stars were observed on 15/2/04 with the exception of HD 183143 and HD 32034 which were observed on 15/9/02 and 
15/8/03 respectively. 



hence phase) at the epoch of observation, and displays the re- 
sults of the spectropolarimetry. These results are discussed in 
Section 4. 

2. Observations & Data Reduction 

The linear spectropolarimetric data were taken during four ob- 
serving runs. Three were done at the 3.9m AAT using the 
25cm camera of the RGO Spectrograph in September 2002, 
February 2003 and August 2003. The fourth was done at the 
4.2m WHT using ISIS in December 2003. The spectroscopic 
data was taken on the AAT runs as part of bad-weather back-up 
programmes, and exactly the same set-up was used as the spec- 
tropolarimetry data, minus the polarization optics. The weather 
was mixed during each run, with some cloud around on some 
of the nights data was taken. The seeing varied from sub- 
arcsecond to > 2" conditions. Additionally, the AAT archive 
was searched for previous spectropolarimetric observations of 
LB Vs. Only high/medium resolution, previously un-published 
data was selected. 

The instrumental set-up was similar for all the AAT spec- 
tropolarimetry observations. A dekker mask with two holes of 
size 2.7" and separation 22.8" was placed in front of the slit 
in order to observe the target and the sky simultaneously. A 
half-wave plate is used to rotate the polarization of the incom- 
ing light, and a calcite block splits the light into two perpen- 
dicularly polarized beams (the O and the E rays), hence four 
spectra are recorded - the O and E beams, of both the target 
and the sky. One complete polarization observation consists of 
one exposure at four different waveplate positions - 0° and 45° 
(to measure Stokes Q), 22.5° and 67.5° (to measure Stokes U). 
Several sets of each object were obtained to check for repeata- 
bility of results, and also to avoid saturation of the CCD at the 
peak of the Ha line. Spectropolarimetric and zero-polarization 
standards were observed each night, and calibration spectra 
were taken after each target by observing a copper-argon lamp. 



The August 2003 and February 2003 runs used a 1200R 
grating in conjunction with the eev2 CCD windowed to 
200x2800, giving a spectral range of 1250A centred at 6500A. 
The September 2002 run used a 600V grating in conjuction 
with the mitll3 CCD windowed to 200x3584, giving a spec- 
tral range of 2690A centred at 5600A. P Cygni was observed 
in December 2003 at the WHT, and the set-up is described in 
Vink et al. (2003). The marconi2 CCD was used in conjunc- 
tion with the 1200R grating, which gave a range of 1055 A cen- 
tred on 6500A. A log of the observations, including integration 
times and spectral resolution achieved is shown in tables 1 and 
2. 

The reduction steps included bias subtraction, flat-fielding, 
cosmic-ray removal, chip linearity correction, spectrum extrac- 
tion and wavelength calibration. The sky spectra were sub- 
tracted from the object spectra in order to compensate for at- 
mospheric features and extended nebulous emission which may 
contaminate the star aperture. These steps were done using the 
figaro package maintained by Starlink. Determination of the 
Stokes parameters was done using the Time-Series/Polarimetry 
(tsp) package, also maintained by Starlink. From the Stokes pa- 
rameters, the degree of polarization P and position-angle (PA) 
9 were found from, 

p = ve 2 + u 2 (i) 

1 u 

6=- arctan — (2) 

PA calibration and wave-plate chromicity were corrected 
for using the Starlink package polmap, by fitting a 4th-degree 
polynomial to the spectra of polarized standards. The high- 
frequency 'ripple' observed in high S/N spectropolarimetry 
data, caused by internal reflections within the waveplate as- 
sembly (Aitken & Hough 2001; Harries & Howarth 1996), was 
corrected for using a Fourier transform algorithm, also within 
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polmap. Instrumental polarization, determined from the zero- 
polarization standards, was found to be less than 0.1%, with 
the origin in Q - U space within the lcr rms scatter. As we 
are only interested in the internal precision of the instrument in 
order to detect changes in polarization across spectral features, 
this has not been corrected for. 

The residuals between the measured Stokes parameters of 
polarized standards and their literature values showed a lcr 
scatter around zero of ~ 0.2%. The uncertainty in position- 
angle Ad ~ 29° x AP/P, and is around 1.5° for P = 4% and 
around 6° forP = 1%. 

There were several extra steps required in the reduction. 
The O and E spectra, as well as spectra from different wave- 
plate positions, were found to be spectrally misaligned with 
each other by a few tenths of a pixel. As the determination 
of the Stokes parameters essentially requires the spectra be di- 
vided by one another, this introduces artifacts in the polariza- 
tion spectrum around sharp features. To compensate for this 
the spectra were cross-correlated with each other and realigned 
according to the measured shift. 

Also, the O and E spectra have slightly different spectral 
resolutions, again by a few tenths of a pixel, despite efforts dur- 
ing focus set-up. We suspect that this is responsible for the er- 
ratic behaviour in the polarization spectrum in unresolved spec- 
tral features, such as a narrow P Cygni absorption components, 
which is highlighted by comparing the different sets of an ob- 
ject. Degrading the spectral resolution to make all spectra the 
same was not found to solve this problem. We suggest this is 
due to pixelation inherent in the data due to the finite separa- 
tion of pixels on the CCD. We identified such artifacts by look- 
ing at the behaviour around narrow features in each data-set, 
and have either edited them out through linear interpolation or 
have been flagged in the text. This problem is unique to studies 
of objects with narrow features such as the discrete absorption 
components in LBV emission lines; it did not affect the results 
of other objects studied in previous papers (e.g. Oudmaijer & 
Drew 1999; Vink et al. 2003). 
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Temperature (K) 

Fig. 1. The equivalent widths as a function of temperature for three 
photospheric absorption lines of early-type supergiants. Circles indi- 
cate stars with emission lines in their spectra, crosses stars with no 
emission. Note that the emission and non-emission stars follow the 
same relation. The uncertainties are of order the size of the plotting 
symbols and are not shown. The temperatures of the spectral types 
were obtained from Schmidt-Kaler (1982). 



3. Results 

The spectra of all the stars show strong Ha emission, often with 
P Cygni absorption components. Many show emission lines of 
Fe ii and [N n], whilst He i and Si n are seen in both emis- 
sion and absorption. Of the objects for which more than one 
set exists, only AG Car exhibits significant spectral variability. 
The spectrum has changed from very reminiscent of HR Car to 
one with emission lines only, indicating a significant change in 
temperature and/or wind properties between the two epochs. 

In Section 3.1 we estimate the temperature, and hence 
phase (whether at minimum or maximum) at the epoch of each 
observation. In Section 3.2 we display the results of the spec- 
tropolarimetry and discuss each object individually. 

3.1. Temperature and Phase at Observation 

The properties of line-driven stellar winds vary with effective 
temperature for two reasons - firstly, as a star cools a growing 
mismatch is formed between its flux-peak and the bulk of the 



driving lines in the UV (Lamers & Cassinelli 1999). Secondly, 
at around 21,000K the opacity dramatically rises as the recom- 
bination of Fe iv to Fe m increases the number of opacity- 
enhancing ions (the so-called 'bi-stabilityjump', see Vinket al. 
1999). There is also evidence for a second bi-stabilityjump at ~ 
10,000K, where Fe m recombines to Fe n (Lamers et al. 1995; 
Vinketal. 1999). 

As LBV effective temperatures can vary from 9,000K to 
up to 30,000K on timescales of a few years, we may expect to 
find different wind properties at different epochs. If the star's 
temperature is latitude-dependent, due to e.g. rotation, we may 
expect to find differential wind properties between equator and 
pole as the star approaches its bi-stability jump, and there- 
fore different wind geometries for the same star depending on 
the phase of the star at observation. This mechanism has been 
used to explain the two-component wind of B[e] supergiants 
by Pelupessy et al. (2000). Here we estimate the temperature 
of the star at observation using two methods - comparing the 
photospheric spectral absorption features with an atlas of early- 



Ben Davies et al.: Asphericity and dumpiness in the winds of Luminous Blue Variables 



5 



Star Log (7W) 





BC method 


Abs. Line method 


AG Car (6/5/94) 


3.98 


±0.03 


4 00 + 04 


AG Car (10/2/03) 


4.25 


±0.03 


_ 2 


HR Car (16/3/92) 


3.87 


±0.05 


3.98 ±°° 2 

0.03 


HR Car (10/2/03) 


3.92 


±0.03 


3.95 

0.05 


7] Car 




_i 




Hen 3-519 




_3 


2 


YX7"D A HZ. 1 

WKA /j1 


4.47 


±0.03 


A /in A £0 

4.4U .. 4.JZ 


PCyg 




_7 


1 


HD 160529 


4.00 


±0.03 


3 95 + 04 

3 - yJ x 0.05 


R40 


4.04 


±0.08 


3 98 + 03 

J - y ° -0.04 


SDor 4 


> 


3.93 


3 95 + - 03 

J - yJ -0.04 


R71 


4.04 


±0.12 


4 19 + 003 

^■ ly -0.04 


R 110 


4.04 ± 0.08 8 


4 00 + 02 


R 116 


4.19 


±0.04 


> 4.15 6 


R 127 


4.29 


±0.04 


413 + 008 
J -o.n 


R 143 






3 99 + 04 

-0.05 



Table 3. Effective temperature at observation as calculated from the 
BC and spectral characteristics methods (see text for details). The two 
estimates agree well for most of the objects, exceptions are explained 
below. 

1 The poor understanding of the intrinsic brightness of 77 Car means 
that the BC method does not produce a reliable result. 

2 The star's spectrum does not contain the diagnostic absorption 
lines. 

3 An insufficient amount of data is available on the star's brightness 
history and brightness at the time of observation. 

4 There is disagreement on the star's temperature at minimum and 
maximum (as stated in vGOl), giving two possible values of the star's 
temperature at observation using the BC method. Also the star's bolo- 
metric luminosity probably changes during outbursts (vGOl), invali- 
dating the BC method. 

5 No AAVSO data exists on this object. 

6 The Si 11 absorption lines appear to be filled in by emission in this 
object, so the temperature was determined from the He 1 line only. 

7 P Cyg is not presently observed to be variable, therefore the BC 
method is not applicable. 

8 R 110 probably does not vary at constant L BO l, and so the BC 
method is invalid (vGOl). 



type supergiants, and via the stars' bolometric correction from 
recent light-curves. 

3.1 .1 . A Spectral Atlas of Early-Type Supergiants 

The spectra of emission line and non-emission line early- 
type supergiants in the range 6240 - 6700A can be found in 
Appendix A, and the full data is available electronically at 
http : //vizier . u-strasbg . fr/viz-bin/VizieR. Among 
the various spectral features, three lines in particular display 
a strong temperature dependence - the Si n AA6347, 6371 
lines and the He 1 /16678 line, whilst no difference is observed 
between the emission and non-emission line stars. Figure 1 
shows the measured equivalent widths (EWs) of these lines 
as a function of effective temperature. The EW of the He 1 
line clearly increases with increasing temperature, appearing at 



around 8000K (spectral-type A7) and peaks at 19 - 21,000K 
(Bl - B1.5), in agreement with Schmidt-Kaler (1982) who 
states that He 1 peaks at spectral type B0. The Si 11 lines de- 
crease in EW with temperature, peaking at 9 - 10,000K (B9 
- A2) and disappearing at around 19,000K (B2). This agrees 
well with Schmidt-Kaler (1982), who states that Si n lines 
peak at spectral-type AO (9700 K). We find no difference be- 
tween emission and non-emission stars. It appears, therefore, 
that these three lines are reasonable indicators of effective stel- 
lar temperature in the range 9,000K - 20,000K (~ A2 - Bl). 
By comparing with the spectral features of the LBVs we can 
obtain an empirical estimate of r e ff, and therefore of the phase 
of the LBV at observation. This empirical relation will break 
down if the lines are in emission, but will only give a mislead- 
ing temperature if there is a tiny amount of emission which 
only partially fills in the absorption. 

3.1.2. The Bolometric Correction Method 

The light variations of LBVs occur at approximately constant 
bolometric luminosity (van Genderen 2001, VGOl hereafter). 
Their apparent variability is caused by their spectral energy 
distribution (SED) shifting in and out of the visual band as 
their effective temperature changes. If we know the apparent 
bolometric magnitude, we can use the difference between this 
and the apparent magnitude at observation (the bolometric cor- 
rection, BC) to determine the star's temperature. This method 
has previously been used to determine LBV temperatures by 
Lamers et al. (1998). 

The BC as a function of r eff for supergiants is quoted in 
Schmidt-Kaler (1982). We have fitted this to within 0.05mags 
with the following relations: 

BC = +27.22 - 6.74 log T eS , (log r eff > 4.25) (3a) 
BC = +15.47 - 3.98 log r eff , (log r eff < 4.25) (3b) 

VGOl quotes the my at visual minimum, with the asso- 
ciated r e jf, for each LBV as collated from multiple sources. 
From this, the BC and hence m^ox can be calculated. The my 
at observation was obtained from the American Association 
of Variable Star Observers (AAVSO, http://www.aavso.org/; 
Waagen 2005), and the BC at the epoch of observation was 
calculated from the difference between this and meoi- Equations 
(3a) and (3b) were then used to calculate the r e jf at observation. 

The estimations of effective temperature at observation are 
shown in Table 3. The two independent measures of tempera- 
ture agree well for each object. Exceptions are marked in the 
table, and probable explanations given in the footnotes. Figure 
2 illustrates the phase of each LBV at observation by placing 
them on a HR diagram, along with their locations at visual min- 
imum and maximum. 

3.2. Wind Geometry 

As explained in Section 1, changes in polarization across emis- 
sion lines supply evidence for aspherical stellar winds. The 
strength of the change in polarization is dependent on both the 
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Fig. 2.HR diagram illustrating the 
temperature of each LBV at the 
epoch of observation. Filled and open 
circles show the temperature deter- 
mined at observation for LBVs and 
candidate LBVs respectively, with 
the solid lines showing the error 
in each measurement. The luminosi- 
ties, and minimum and maximum 
observed temperatures were collated 
from the literature by Smith et al. 
(2004). The Humphreys-Davidson 
limit is marked by a dashed line. 
r\ Car, Hen 3-519 and P Cyg are 
marked on the plot despite the fact 
that their temperatures could not be 
determined - values are taken from 
Smith et al. and Cox et al. (1995). 



density contrast between different regions in the wind and the 
inclination angle. A non-detection could imply either the wind 
is symmetric about some axis (e.g. an equatorially-enhanced 
wind) but is oriented face-on; or the electron density at the base 
of the wind is isotropic (i.e. the wind is spherically symmet- 
ric) to within the detection limit. This detection limit is depen- 
dent inversely on the S/N of the spectrum, and the contrast of 
the emission line as the line-emission is 'diluting' the polar- 
ized flux from the continuum. The detection limit per pixel for 
the maximum intrinsic polarization AP intMax for each object is 
given by: 



AP mt , max (%) = J— x j^- x 100 (4) 



where L is the line-to-continuum contrast. In the cases where 
no line-effect is observed, we state the systematic detection 
limit of our observations. Where line-effects are detected, the 
PA of the intrinsic component of the star's observed polariza- 
tion is determined from the line-to-continuum vector in Q - U 
space using Eqns. (1) and (2). Note that this is independent of 
the interstellar polarization (see below). 

The continuum polarization is discussed below, and proper- 
ties of the Ha emission as well as continuum polarization mea- 
surements are shown in Table 4. Object-by-object descriptions 
are listed below, starting with the Galactic objects. Of the 14 
targets observed, 6 show definite line-effects, with one further 
possible detection. Of the 4 objects for which more than one 
observation exist, 3 show line-effects, and each of these three 
display polarimetric variability. 



3.2.1. Continuum Polarization 

The observed polarization in the continuum is a superposition 
of multiple sources. We assume that this is broadly due to three 
mechanisms: scattering of light by free electrons at the base of 
the stellar wind; scattering by gas and dust in the star's neb- 
ula; and dichroic absorption by dust in the interstellar medium 
(ISM). Here we focus on the polarizing effects of asphericities 
in the stellar wind within a couple of stellar radii, which can 
vary on short timescales (days). For the purposes of this paper 
we will lump the nebula and the ISM together, and say the con- 
tinuum polarization is due to the interstellar polarization (ISP) 
and electron scattering in the stellar wind. 

Determination of the ISP, and hence separation of the two 
polarizing components is not straight-forward. The field-star 
method, attempted by Schulte-Ladbeck et al. (1994); Harries 
et al. (1998); Parthasarathy et al. (2000) amongst others, gen- 
erally yields mixed results due to uncertain distances, the an- 
gular scale on which the ISP can change, and the unknown in- 
trinsic polarization of the field stars. The line-centre method, 
used by Schulte-Ladbeck et al. (1993, 1994); Clampin et al. 
(1995); Oudmaijer et al. (1998), assumes that as the Ha line- 
emission is intrinsically unpolarized and formed over a much 
larger volume it will undergo negligible scattering and re- 
main unpolarized. Therefore, if the line is much stronger than 
the continuum (such that the line contrast term in Eqn. (4) 
tends to 1), the Ha line-centre polarization is a measure of the 
ISP. Multiple Spectropolarimetric measurements of AG Car by 
Schulte-Ladbeck et al. (1994) showed that the polarization of 
Ha remained constant over ~ 1.5 years whilst the continuum 
polarization varied, showing that this method can be an effec- 
tive method of determining the ISP. 
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Table 4. The observed Ha data. The first column denotes whether or not a line-effect was detected; the second column shows the contrast 
between the line-peak and the continuum to an accuracy of 0.1; the equivalent-widths and FWHM are measured to an accuracy of ~ 5%. 
The continuum polarization was measured from featureless sections either side of the emission line. As the R 40 data from consecutive nights 
showed no variability, the data was combined to improve the S/N. The quoted uncertainties do not include the external errors, which we estimate 
to be 0.1 - 0.2% (see Section 2). 



3.2.2. Object-by-Object Descriptions 

AG Car One of the prototype LBVs, AG Car has been highly 
variable for over 30 years with light variations of over 2 mags. 
It is surrounded by an elliptical nebula which has a hint of a bi- 
polar morphology, and a dynamical age of ~ 10 4 years (Nota 
et al. 1992). High-resolution HST images have shown that the 
nebula is made up of individual clumps, filaments and bubbles 
(Nota et al. 1995). Multiple spectropolarimetric observations 
revealed evidence for an axisymmetric stellar wind, with a PA 
alternating between perpendicular and parallel to the major- 
axis of the nebula (Schulte-Ladbeck et al. 1994, SL94 here- 
after). SL94 attributed this to either a wind which flips between 
an equatorial- and polar-enhanced wind, or multiple-scattering 
effects due to changes in optical depth. 

The data from both epochs show line-effects (see Fig. 3). 
The archive data from 1994 had problems with varying spec- 
tral resolutions from set to set. Consequently, erratic polariza- 
tion behaviour around the emission line was observed, and has 
been edited out in order to show the broad PA rotation coin- 
cident with the broad wings. Similar behaviour is seen in the 
2003 data, but the continuum PA has changed by ~ 30°. The 
varying PA of the observed polarization alone is evidence for an 
intrinsic component, as the interstellar polarization (ISP) will 
not vary on these timescales (see Sect. 3.2.1). 



The line-centre polarization of AG Car has previously been 
shown to remain roughly constant, and is therefore probably a 
good estimate of the ISP (SL94, see Sect. 3.2.1). As shown 
in Fig. 3, the ISP measured by SL94 agrees well with the 
line-centre polarization of our 2003 data. The full Q - U ex- 
cursion of the 1994 data has been attenuated due to the line- 
centre being edited out, but it can clearly be seen to point 
towards the same region of Q - U space. The mean of the 
line-centre polarizations, including those measured by SL94, 
is Q = 0.36 + 0.07; U = -1.08 + 0.09. 

Due to the attenuation of the 1994 data, the intrinsic po- 
larization is measured from the line-wings, and is found to be 
P > 0.26 + 0.03%, 6 = 26 + 3°. The intrinsic polarization of the 
2003 data is P = 0.41 ± 0.03%, 6 = 84 + 2°. A weak line-effect 
is also seen across the He i /16678 line which is in emission at 
this epoch, with a similar PA to Ha. However, line-effects are 
not observed across any other emission lines. The PAs from 
both epochs, 26° and 84°, appear to be aligned with neither the 
nebula (~ 135°), nor the jet (~ 35°). 

Looking at our data and the measurements of SL94, it can 
be seen that the intrinsic polarization of AG Car varies greatly 
from epoch to epoch (Fig. 3). The data from SL94 showed the 
intrinsic polarization vector flipping from one side of the ori- 
gin in Q - U space to the other, leading them to suggest that 
the intrinsic polarization varied between two preferred planes 



8 



Ben Davies et al.: Asphericity and dumpiness in the winds of Luminous Blue Variables 



AG Car 6/5/94 err=0.05% 



AG Car 10/2/03 err=0.07% 




6500 6550 6600 6650 
Wavelength (A) 




Polarization Measurements of AG Car 



15/3/92 





6/5/94 



10/2/03 



) Ha line-centres 



(£) 12/2/93 
'fy 27/11/91 



0.5 
Q/I (%) 



6500 6550 6600 6650 
Wavelength (A) 



Fig. 3. Left, Centre: Polarization spectra of AG Car from the two dates observed. The bottom panel of each triplot shows the Stokes /, or 
intensity spectrum; and the middle and top panels show the degree and PA of the polarization respectively as a function of wavelength. The 
data has been rebinned to the lcr error indicated above each plot. The 6/5/94 data displayed erratic behaviour around the absorption component 
and line-centre due to reasons explained in Sect. 2, and was edited out of this plot for clarity. 

Right: Multiple spectropolarimetric measurements of AG Car plotted in Q - U space. Circles are continuum points from SL94; the cross is 
from archival data; and the square is our continuum data. Dotted lines show the excursions to the Ha line-centre polarization. Note that the 
1994 excursion was attenuated due to fluctuating data (see Section on AG Car), and that SL94 averaged their three line-centre polarizations 
together. The mean of the line-centre polarizations is marked on the plot, with a lcr rms error of the size of the symbol. Uncertainties in the 
continuum are limited by the external errors, of order 0.1%. 



90° to each other. These planes aligned with the major and mi- 
nor axes of the coronographic image of the AG Car nebula in 
Clampin et al. (1993). Our observations lie in neither of these 
planes, with the archive data being taken only ~ 6 months af- 
ter the last SL94 observation. The continuum measurements of 
Leithereret al. (1994, not shown in Fig. 3) are different again. If 
the wind was indeed flipping from an equatorial to a polar flow, 
we would expect the polarization of any intermediate phases to 
lie along the same axis in Q - U space but with a smaller line- 
to-continuum vector. This is not observed here, with strong de- 
polarizations observed at all epochs at unrelated PAs. Unless 
the assumption used to estimate the ISP are invalid, the evi- 
dence points away from the flip-flopping axi-symmetric wind 
scenario. Instead, it may point towards a 'clumpy' wind model, 
where the PA of the intrinsic polarization at any time is re- 
lated to the projected angle between the star and the dominant 
clump(s). This is the same scenario to that proposed for P Cyg 
by Nordsieck et al. (2001, see section on P Cyg). 

H R Car Like AG Car, HR Car has been strongly variable for 
over 25 years (> 1 mag on timescales of ~ 3 years). It is sur- 
rounded by a filamentary bi-polar nebula (PA ~ 135°) which 
has a kinematic age of ~ 5000 years (Nota et al. 1997). The star 
was noted by Serkowski et al. (1975) as polarimetrically vari- 
able, implying it has intrinsic polarization. This was confirmed 
by subsequent polarimetric monitoring by Parthasarathy et al. 
(2000). It was observed spectropolarimetrically by Clampin 



et al. (1995), who found evidence for intrinsic polarization at 
a PA of 30 ± 1°, which is aligned with the 'finger' of nebu- 
lar emission (PA ~ 35°) within ~ 3" of the central star noted 
in their coronagraphic image. They also noted the presence of 
a number of clumps of emission close to the star, which did 
not correspond to the filamentary structures of the nebula and 
appeared to be roughly spherically distributed around the star. 

In our spectropolarimetry, both epochs show line effects 
(see Fig. 4). The 1992 data shows a depolarization accom- 
panied by a PA rotation. The 1992 intrinsic polarization is 
P = 0.67 + 0.03%, 6 = 138 ± 2°, which is aligned parallel 
to the major axis of the bi-polar nebula. The 2003 data shows 
a lower continuum polarization, with an increase in polariza- 
tion across the emission line. The 2003 intrinsic polarization 
is P = 0.56 ± 0.03%, 0=13 + 2°. This PA bears no obvious 
relation to the PA measured at the previous epoch, nor to any 
large-scale nebular features. 

The line-centre polarization from our observation, the 
archival data and the data in Clampin et al. (1995) agree well 
with each other, indicating that this is a good estimate of the 
ISP towards HR Car at 6563A (as with AG Car). The mean of 
these measurements is Q = -0.80±0.05%, U = -3.56±0.03%; 
implying P = 3.66 + 0.05%, 6 = 128.7 + 0.4° (see Figure 4, 
right). 

Line-effects are also seen across Fe n lines at both epochs. 
In the 1992 data, there is a borderline detection across the 
i6456 line, with a line-to-continuum PA of ~ 140°, i.e. par- 
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Fig. 4. Left, Centre: Polarization spectra of HR Car from the two dates observed. The panels show the same as Figure 3, and the data has been 
rebinned to a la error of 0.05%. 

Right: Multiple spectropolarimetric measurements of HR Car, plotted in Q - U space. The cross marks the continuum polarization at Ho- for 
the 1992 archival data; the circle for the data taken by Clampin et al. (1995); and the square for our data. The dotted lines show the excursions 
to the line-centre polarization. The mean of the line-centre polarizations is marked, with a la rms error of the size of the symbol. Uncertainties 
in the continuum are limited by the external errors, of order 0.1 - 0.2%. 



allel to the Ha depolarization. In the 2003 data, line-effects are 
detected across the Fe n (AA 6417, 6433, 6456, 6516) lines. 
The depolarizations all have PAs between 35 - 60°, but these 
are borderline detections. 

Taking the ISP to be that calculated above we find that the 
Clampin data and archival data, taken approximately 1 year 
apart, are perpendicular to each other with PAs of ~ 45° and 
~ 135° respectively. These align with the major and minor axes 
of the bipolar nebula. Our data, taken 10 years later, has a com- 
pletely different PA of ~ 15°. Taking the spectropolarimetry of 
Clampin et al. (1995) and multiple /?-band polarimetric mea- 
surements of Parthasarathy et al. (2000) into account, the in- 
trinsic polarization of HR Car seems to have no preferred axis 
with time. In a clumpy wind we may expect this kind of tem- 
poral variability, where the polarization at each epoch is related 
to the projected PA between the star and the closest, densest 
clumps. 

t] Car The most famous of the LBVs underwent a huge erup- 
tion in the mid- 19th century when it increased in brightness by 
over 2 mags and ejected ~ 10M G (Morris et al. 1999; Smith 
et al. 2003b), the product of which is the 'homunculus' neb- 
ula we see today (for an extensive review of this object see 
Davidson & Humphreys 1997). Aside from a smaller eruption 
in the late 19th century, which is possibly responsible for the 
mini-homunculus within the main nebula detected by Ishibashi 
et al. (2003), the star has been more or less stable. The periodic 
(~ 2000 days) spectroscopic events are attributed to the star be- 
ing in a highly eccentric binary system (Damineli et al. 2000); 



and the gradual brightening in the optical and NIR is attributed 
to the expansion of the nebula (Davidson & Humphreys 1997). 
t] Car's wind as a function of latitude was studied by Smith 
et al. (2003a) by looking at starlight reflected of the surround- 
ing nebula. They determined that the wind has higher terminal 
velocities and higher mass-fluxes towards the poles, assuming 
the star's axis is aligned with the bipolar nebula. Inteferometric 
observations by van Boekel et al. (2003) revealed that the star's 
NIR emission is prolate, with a major axis of PA 134 + 7°. This 
is aligned with the PA of the humunculous nebula (PA 132°, 
Morse et al. 1998). 

Our data shows a very complex polarization profile, corre- 
sponding to an emission line which is made up of three com- 
ponents - the broad wings, the narrower component containing 
the bulk of the flux, and the narrow peak (see Fig. 5). The broad 
wings correspond to the Q - U excursion labelled 1 with a PA 
of 91°; the flux-bulk to the excursion labelled 2 with a PA of 
52°; and the narrow peak to the Q- U loop marked 3, with a PA 
of 95°. There is also a loop as the vectors return from the high- 
flux region to the continuum region. This does not correspond 
to any detectable spectral features. 

Comparison of our spectrum with archive data taken by 
the Hubble Space Telescope (HST) one day previously reveals 
that the narrow emission peak (#3) does not originate from the 
star, but from a nebulous region approximately 1" away. Our 
data does not spatially resolve the star and the nebula, and 
so this nebulous emission will not have been removed at the 
sky-subtraction stage of data reduction (see Sect. 2). The [N n] 
emission line redward of the Ha line, which displays a depolar- 
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ization, also originates from the nebula. As the nebula is located 
well outside the polarizing region close to the star, any polar- 
ization it may have will be unrelated to that of the continuum. 
Closer inspection of the HST spectrum reveals that the nebula 
also has continuum emission, probably from reflected starlight. 

Given that the light we see in our spectrum is made up of 
a direct starlight component and a reflected nebula component, 
these may be responsible for the distinct Q-U excursions. The 
emission component containing the bulk of the flux ('2' on Fig. 

5, which probably corresponds to a line-forming region close 
to the star) with a PA of 52 + 3° is perpendicular to the major- 
axis of the resolved prolate structure detected by van Boekel 
et al. (2003). This is consistent with polarization due to scatter- 
ing off a bipolar wind, as per the scenario suggested by Smith 
et al. (2003a). The level of depolarization, at ~ 1%, is also con- 
sistent with such a scenario (Cassinelli et al. 1987). The other 
Q-U excursions, which seem to have PAs of ~ 95° and a 
much stronger depolarization of > 2%, may be due to scattered 
light off the unresolved nebulous component. This hypothesis 
is supported by other sharp features of the spectrum (not shown 
here), which also show depolarizations with similar PAs of ~ 
90°. These PAs do not relate to any nebular features nor to the 
PA of the slit, but as shown in Schulte-Ladbeck et al. (1999) 
the polarization of 77 Car varies greatly within 2" of the central 
source. 

In summary, the complex polarization profile may be a 
combination of spatially unresolved emission from the nearby 
reflection nebula and emission from the stellar wind. The polar- 
ization of the light from the stellar wind is perpendicular to the 
bipolar axis of the homunculus nebula. This is consistent with 
single scattering off an optically thin bipolar wind, or multi- 
ple scattering off an optically thin equatorially-enhanced wind 
(Angel 1969; Wood et al. 1996). 

Hen 3-519 This object is classified as a post-LBV/pre-WR 
star by (Davidson et al. 1993). It is not observed to be vari- 
able and has broad emission-lines in its spectrum. For these 
reasons it is not possible to determine the star's temperature at 
observation using the two methods outlined in Sect. 3.1. The 
star is surrounded by a roughly circular nebula of diameter ~ 
55". The nebula appears to be a clumpy bubble, rather than a 
detached ring as seen around AG Car (Stahl 1987; Duncan & 
White 2002). The spectra from the two different epochs sepa- 
rated by ~ 9 years look very similar, with the exception of a P 
Cygni profile present in the Ha line of the Feb 2002 data not 
present in the 1994 data. 

The archive data from 6/5/94 (Figure 6, left) does not show 
an obvious line-effect. The data from 10/2/03 shows different 
polarization behaviour across the features of the line (see Fig. 

6, centre). The line-centre polarization is consistent with that of 
the archival data, showing up in Q- U space as a dark cluster of 
points above the continuum region. This behaviour is also seen 
in the Hel AA 5876, 6678 lines. The mean line-to-continuum 
vector, and hence intrinsic polarization, is Q = -0.10 + 0.03%, 
U = -0.14 + 0.03%; P = 0.17 + 0.04%, 6 = 117 + 12°. The 
mean of the line-centre polarizations is P — 2.89 + 0.02%, 8 = 
97.9 + 0.3°, which agrees with the data of the previous epoch. 



Eta Car err = 0.1% 




-5 -4 -3 -2 

q (%) 



Fig. 5. Top: The polarization spectrum of 77 Car. The panels show the 
same as Figure 3. The data has been rebinned to a lcr error of 0.1%. 
The dotted line shows a magnification of the Stokes / spectrum to 
highlight the broad wings and weaker spectral lines. 
Bottom: The spectrum plotted in Q - U space. The arrows join points 
of increasing wavelength. The location of the continuum is marked 
'cont.'. The distinct excursions mentioned in the text are indicated 
with dotted lines. The excursion associated with the wings is marked 
' 1' ; the excursion associated with the narrower region containing the 
bulk of the line-flux is marked '2'; and the loop associated with the 
narrow, unresolved nebular emission is marked '3'. 



As well as the Ha line having a slightly different polariza- 
tion to the continuum, there is a large change in polarization 
in the blue wing of the 2003 data, showing up as a looped ex- 
cursion in Q - U space (see Fig. 6, right). This behaviour in 
the blue-wing is also seen in the He 1 AA 5876, 6678 lines. The 
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Fig. 6. Left, Centre: The polarization spectra of Hen 3-519 from 6/5/94 and 10/2/03 repectively. The panels show the same as Figure 3 and the 
data has been rebinned to the lcr error indicated above each plot. 

Right: The polarization spectrum of the 10/2/03 data plotted in Q - U space. Arrows join points of increasing wavelength. The distinctly 
separate continuum and line-centre regions have been marked, as well as the loop corresponding to the depolarization seen in the blue wing. 



mean change in polarization from the line-centre polarization 
of these loops is P = 0.87 + 0.08%, 6 = 25 ± 3°. At present, we 
have no explanation for these features, which cannot be due to 
the simple dilution of the polarized continuum photons. They 
are also unlikely to be related to the "McLean" effect (McLean 
1979), where an increase of polarization due to scattered light 
is seen across blueshifted P Cygni absorption (see Vink et al. 
2002), as here the features in Hen 3-519 extend beyond the 
blueshifted part of the P Cygni absorption. 

WRA 751 This star was first considered to be an LBV by Hu 
et al. (1990), but is still referred to as an LBV candidate by 
many authors. It is classified as a strong-amplitude variable by 
VG01, but the AAVSO light-curve shows that it has not varied 
by more than 0.5mags over the last ~10 years. Radio images 
by Duncan & White (2002) show a roughly circular, clumpy 
nebula of diameter ~ 20"with evidence for a central, edge-on 
nebular ring of diameter ~ 5". The object's polarization spec- 
trum shows no line-effect (see Figure 7). 

If the light from the star is intrinsically polarized, the upper 
limit we can put on this is rather high due to the low S/N of the 
spectrum and the relative weakness of the Ha line (see Sect. 
3.2, Eq. (4) ). The sensitivity limit for detection of dilution of 
the continuum by the emission line is 0.44%. As we have pos- 
itive detections in this study of less than 0.3%, we cannot rule 
out a similar level of asphericity in WRA 75 1 's wind. 

P Cyg Since undergoing two huge eruptions in the 17th cen- 
tury, there is no evidence this star has undergone any short- 
timescale variability, although de Groot & Lamers (1992) ar- 
gue the star has gradually brightened by about 0.5 mags since 
then as it evolves redward. The star's nebula has been the sub- 



ject of much investigation over the last ~ 10 years. It consists 
of an inner- and outer-shell, of diameters 22" and 1 .6' respec- 
tively, both with roughly spherical morphologies (Barlow et al. 
1994; Meaburn et al. 1996). Optical images of the circum- 
stellar material reveal a complex clumpy morphology down 
to the inner 0.2" (Meaburn et al. 2000; Chesneau et al. 2000). 
Multiple interferometric radio observations have shown that the 
dumpiness of the inner 0.4" of the nebula is highly variable on 
timescales of days (Exter et al. 2002). 

Both the observations of this object, taken 3 days apart, 
show changes in polarization just blueward of the centre of the 
emission line, with varying behaviour from set to set. They line 
up not with the P Cygni absorption component but with the 
steep rise in flux from the absorption component to the centre of 
the emission line, which may be at the limit of the instrument's 
spectral resolution. As polarization in this part of the spectrum 
varies wildly from set to set, whilst the continuum level remains 
constant, we therefore explain these features as artifacts intro- 
duced by the slightly different spectral resolutions of the O and 
E beams emerging from the calcite. In terms of difference in 
polarization between the continuum and the line-centre, there 
is no change within the uncertainties (see Fig. 7). The detection 
limit is 0.15% and 0.18% for the 10/12/03 and 13/12/03 obser- 
vations respectively. The Ha line-centre polarization we mea- 
sured agrees well with the line-centre polarization measured by 
Taylor et al. (1991) and Nordsieck et al. (2001), supporting the 
hypothesis that the Ha line is intrinsically unpolarized. We find 
no change in continuum polarization for our two datasets taken 
3 days apart. 

Multiple spectropolarimetric measurements of P Cyg are 
described by both Taylor et al. (20 measurements) and 
Nordsieck et al. (15 measurements). Both found that while the 
line-centre polarization remained constant the continuum po- 
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Fig. 7.Polarization spectra of the Galactic 
LBVs WRA 751, P Cyg and HD 160529; 
and the SMC LBV R 40. The panels show 
the same as figure 5, and the data has been 
rebinned to the lcr errors stated above each 
plot. Note the different wavelength range of 
P Cyg. The jump in polarization seen in the 
Ha absorption component of P Cyg is at- 
tributed to the steep increase in Stokes / be- 
ing unresolved, as explained in section 2. 
HD 160529 did not display a line-effect at 
either epoch it was observed, and only the 
Feb 2003 data is shown here. 



larization was variable on timescales of days. Both interpreted 
this as variable intrinsic polarization of the object coupled with 
the ISP, with an intrinsically unpolarized Ha line as described 
above for AG Car and HR Car. Taylor et al. found that the in- 
trinsic polarization varied from 0.04% to 0.48% at random PAs. 
Such behaviour was ascribed by Nordsieck et al. to electron 
scattering by a clumpy wind, with a spherical distribution of 
clumps and a clump-ambient material density contrast of > 20. 
It may be the case that at the epoch of our observations any 
intrinsic polarization due to wind clumps was too small to be 
detected. 

HD 160529 Despite being established as spectroscopically 
and photometrically variable for over 30 years (Wolf et al. 
1974), Serkowski et al. (1975) noted no polarimetric variabil- 
ity, and it is listed by Hsu & Breger (1982) as a 'polarized stan- 
dard' . This star represents the lower-luminosity end of the LBV 
instability strip and is not strongly variable, with photometric 
variations of only a few tenths of a magnitude over the last 20 
years (Sterken et al. 1991 ; Stahl et al. 2003). The star has no as- 
sociated nebulosity, no nebulous emission in its spectrum and 
no evidence of circumstellar dust (Hutsemekers 1997; Nota & 
Clampin 1997). We note no [N n] emission in our spectra. 



The relative weakness of the Ha emission makes the detec- 
tion of a line effect difficult for this object (Fig. 7. The detection 
limit for a line-effect for this observation is 0.4%. 

R 40 This star is the only known LBV in the SMC. Sterken 
et al. (1998) show that R 40 increased steadily in brightness by 
0.7 mags in the years 1985 - 1995. AAVSO records show this 
object has been constant in luminosity since ~ 1995, although 
there is a gap in their records of ~ 3 years (Waagen 2005). 
There is presently no evidence for associated nebular or dust 
emission, and we find no [N n] emission in its spectrum. 

No line-effect was observed on either of the observing 
dates. As the object was observed on consecutive nights and 
the data for each night looked the same, the data were com- 
bined to improve the S/N, but still no line-effect was observed 
(Fig. 7. The detection limit is 0.37%. 

S Dor Despite being one of the 'archetypal' LBVs for many 
years, the star has only recently been observed with the spec- 
tral appearance of an F-type supergiant - the common spectral- 
type of LBVs at visual maximum (Massey 2000). Also, it has 
no directly detectable LBV nebula. Weis (2003) found that the 
dynamical age of the large-scale elliptical nebula surrounding 
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Fig. 8. Polarization spectra of the LMC LBVs S Dor, R 71, R 110, R 116, R 127 and R 143. The panels show the same as figure grrrr. Note 
the different wavelength scale for R 127. Where line-effects are observed, the Stokes / spectrum is magnified to show the broad change in 
polarization associated with the wings of the emission line. A narrow P Cygni absorption component in the Ha line of R 127 was found to 
introduce artifacts in the polarization spectrum (see Section 2), and has been edited out. 



S Dor meant that it was too old to have been formed in the 
LBV phase; and despite broadened [N n] emission in the star's 
spectrum (also noted in our spectrum), any circumstellar nebula 
from which it originates would have to be less than the detec- 
tion limit of ~ 0.25pc in diameter. The star's spectrum indicates 
a T e ff at observation of ~ 9000K, putting the star close to visual 
maximum. The bolometric correction (BC) method of deter- 
mining the star's temperature is unreliable, due to disagreement 
over the star's temperature and luminosity at visual minimum, 
and also that the star may vary in bolometric luminosity during 
outbursts (see VG01). The AAVSO light-curve of S Dor shows 
it to have been exhibiting strong-amplitude variability over the 
last 10-15 years. 

The polarization spectrum (Fig. 8) has a suggestion of a 
broad PA rotation across the line from ~ 60° to ~ 40°. This cor- 



responds to a 2cr detection and is therefore a borderline case. 
The single pixel dip in PA at the centre of the line is not trusted 
as 'real' behaviour, as sharp spectral features can produce er- 
ratic data for the reasons described in Sect. 2. The detection 
limit in this data is 0.2%. 



R 71 The star is noted by Weis (2003) as having no visually 
detectable nebula and no nebular emission in its spectrum, al- 
though Voors et al. (1999) did find evidence for a warm dust 
component. Also, our spectrum does show very weak [N n] 
i6584 emission which suggests a nebular component. VG01 
states that this star has never been observed as cool as the clas- 
sical LBV minimum temperature of ~ 8500K. AAVSO records 
show this star has been constant in brightness to within ~ 0.2 
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mags for the last ~ 15 years, although the star did increase in 
brightness by 0.5mags for about a year in 2003 (Waagen 2005). 

This object shows no line-effect (Fig. 8). The comparitive 
weakness of the Ha line in this data means that the detection 
limit is 0.4%. 

R 1 1 It is not known whether this star has associated circum- 
stellar dust or gas as it is poorly-studied. We do find, however, 
weak [N n] /I6584 emission in our spectrum, which suggests 
that the star may have some form of nebula. The star was ob- 
served to increase in brightness by 0.9 mags to m v = 9.7 be- 
tween 1960 and 1993 by van Genderen et al. (1997). AAVSO 
photometric data on this object is limited, but it has been seen 
to decrease in brightness to m v ~ 11 in the last ~ 10 years 
(Waagen 2005). 

The data exhibits a broad depolarization of > 0.5% accom- 
panied by a broad PA rotation (Fig. 8). The width of this fea- 
ture corresponds to the width of the broad wings of the emis- 
sion line. The sharp features observed in the polarization spec- 
trum at the line centre are spread over ~ 10 pixels, which imply 
that they are resolved and are judged to be real. The line-to- 
continuum vector in Q— U space yields an intrinsic polarization 
of (0.55 + 0.09%) at a PA of (26 ±5)°. 

R 1 1 6 The similarities between this star and £ i Sco noted by 
van Genderen & Sterken (1999) led them to suggest that R 1 16 
is an ex/dormant-LBV. AAVSO records show that the star has 
been constant in brightness over the last 15 years to within 0.2 
mags (Waagen 2005), and photometry collated from various 
sources by van Genderen & Sterken suggested the star has been 
constant for over 40 years. There is presently no evidence for 
associated nebular or dust emission. 

This object shows no line-effect (Fig. 8). The detection 
limit is 0.28%. 

R 127 The LBV nature of this star was discovered by Stahl 
et al. (1983). Since then it has been observed to cool to the clas- 
sical LBV minimum temperature of ~ 8000K around 1993, and 
is now apparently on its way back to visual minimum (VG01 
& references therein; Waagen 2005). Coronagraphic imaging 
by Clampin et al. (1993) revealed an axisymmetric circumstel- 
lar nebula, with density enhancements either side of the star 
at a PA of ~ 100° and an opening angle of 90°. HST observa- 
tions by Weis (2003) reveal evidence for asymmetric expansion 
within the inner nebula. 

Spectropolarimetry by Schulte-Ladbeck et al. (1993) re- 
vealed a high degree (~ 1%) of intrinsic, variable polarization; 
the PA of which remained roughly perpendicular to the density 
enhancements. They explained the roughly constant intrinsic 
PA and observed P Cygni profiles as electron-scattering off an 
expanding equatorially-enhanced wind aligned with the outer 
nebula. The variable degree of polarization and slightly vary- 
ing PA were explained by a dumpiness of the wind. 

Our data shows a broad rotation in PA which corresponds 
to the broad wings of the emission line (Fig. 8). The emission 
line has a narrow P Cygni absorption component which was 



found to introduce erratic behaviour in the polarization spec- 
trum. For the sake of clarity this feature was edited out. The 
sharp increase in PA at the line-centre is spread over ~ 5 pixels 
and so are judged to be real. The line-centre polarization has 
been attenuated slightly by editing out the absorption compo- 
nent, but the Q - U excursion approaches the Ha line-centre 
polarization measured by Schulte-Ladbeck et al. (1993). 

Taking this to be the ISP towards R 127, the intrinsic po- 
larization is 0.50 + 0.04% at a PA of 20 + 2°. This is in 
good agreement with Schulte-Ladbeck et al. who found a PA 
of 24.4° in the V band, although at a higher degree of polar- 
ization of 1.37%. These results, along with other polarization 
measurements collated by Schulte-Ladbeck et al., show that the 
intrinsic polarization of R 127 appears to be at a roughly con- 
sistent PA of 25 + 10°. This endorses their proposed model of 
an equatorially-enhanced, clumpy wind viewed at a high incli- 
nation angle. 

R 1 43 Weis (2003) reported that the star is located in a small, 
triangular-shaped nebula of linear size ~ 1.2pc within the 30 
Dor region. The disrupted shape of the nebula was speculated 
by Weis to be due to strong stellar winds, high density and 
turbulent motions present in 30 Dor. When observed by Nota 
et al. (1995), the nebula was found to have filamentary struc- 
tures eminating from the central source, reminiscent of HR Car. 
However, Smith et al. (1998) found that the dynamics of these 
structures suggested they were not associated with the star but 
of the background 30 Dor H n region. They did find evidence 
for a small clump of ejected material within 2" of the star but 
could not draw any conclusions as to its morphology. 

This object shows no line-effect (Fig. 8). The detection 
limit is 0.34% due to poor S/N. 

4. Discussion 

Of the 1 1 confirmed LBVs, we have 5 positive detections of 
line-effects plus the borderline case of S Dor. Of the three LBV 
candidates, only Hen 3-519 shows a change in polarization 
across its emission lines, and the behaviour of this star's po- 
larization spectrum cannot simply be explained by the dilution 
of the polarized flux by the emission line. Of the non-detection 
objects, in some cases the detection limits exceed the levels 
of depolarization observed in objects with line-effects. There 
may therefore be similar levels of depolarization in some of 
these objects that have gone undetected. Even if all LBVs had 
aspherical winds, we would expect some non-detections due 
to inclination and cancellation effects e.g. face-on disk, or two 
clumps at 90° to each other. We therefore put the rate of wind- 
asphericity in LBVs at > 50%, higher than that of the other 
classes of evolved massive star that have been studied compre- 
hensively, O supergiants and Wolf-Rayet stars (25% and 20% 
of sample sizes 20 and 16 respectively, Harries et al. 2002, 
1998). 

Rotation and mass loss play intricate roles in the evolution 
of massive stars (Langer 1997; Maeder & Meynet 2000). Given 
the transitional role of LBVs within the evolutionary paths of 
massive stars, the high incidence of intrinsic polarizations that 
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we find in comparison to those of O/WR stars is expected to 
reveal telling information about the onset of wind asymmetries. 
To probe the origin of the wind asymmetries, we correlate our 
detections with a range of LBV parameters (Sect. 4.1), before 
we discuss possible scenarios in Sect. 4.2. 

4.1. Correlations 

In terms of temperature and phase, no correlation is observed 
between detections and non-detections - line-effects are ob- 
served at visual minimum (HR Car), maximum (R 110), and 
intermediate (e.g. R 127) phases. No correlation is observed 
with luminosity, although we point out that the luminosities of 
these objects (even those in the Clouds) are uncertain, with dif- 
ferent authors stating different values for L/ L . VG01 states 
that luminosities of Galactic objects are uncertain to 0.2-0.3 
dex, while MC objects are uncertain to ~ 0. 1 dex. 

4.1 .1 . Line-strength - wind asphericity relation? 

Inspection of Table 4 and Figure 9 shows that the objects show- 
ing line-effects all have strong Ha lines. Indeed, the strongest 
Ha line definitely showing no line-effect is R 143, with a line- 
to-continuum contrast of 6.4. S Dor has a line contrast of 8.3, 
and has a borderline case for a line-effect. Only P Cyg of the 
strong-emission objects has no clear-cut line-effect. From this 
we could draw one of two conclusions: either the method is 
not sensitive enough to detect line-effects in weak emission 
lines; or objects with strong Ha lines have strong asphericity 
in their winds. Line-effects (albeit weak) are seen across the 
He i /16678 line in the spectrum of AG Car (contrast ~ 1.7), the 
weaker lines in 's spectrum (~ 1.7), and the weak Fe n lines 
in the spectrum of HR Car (~ 1.5). This therefore leads us to 
discount the first explanation and suggest that LBV winds pro- 
ducing stronger emission are more likely to have an aspherical 
geometry. 
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Fig. 9. Strength of depolarization against the line-to-continuum con- 
trast of the Ha line. Circles show line-effects, crosses show non- 
detections. Where no depolarization was observed, the upper limit for 
an undetected line-effect is shown. For the sake of clarity, t] Car is not 
shown as its depolarization and line contrast are both much greater 
than any other object. 
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4.1.2. Galactic vs. Extra-Galactic LBVs 

Of the Galactic objects observed, 4 out of 7 show line-effects. 
Of the three non-detections, one (P Cyg) has been shown to ex- 
hibit a change in polarization across its emission lines at earlier 
epochs (Taylor et al. 1991; Nordsieck et al. 2001). By com- 
parison, only 2 out of 7 Magellanic Cloud objects show line- 
effects, with one borderline case (S Dor). One explanation for 
this difference may be that the selection of Galactic LBVs is 
not truly representative of the LBV phenomenon. Hen 3-519 is 
commonly thought-of as a post-LBV object, indeed its broad 
emission lines are reminiscent of a Wolf-Rayet star. WRA 75 1 
is located in a distinctly separate location on the H-R diagram 
to the other LBVs (see Figure 2) and has not been observed 
to show variability on the same scales as other LBVs such as 
AG Car or S Dor (van Genderen et al. 1992). r] Car and P Cyg 
are unusual objects even by LBV standards - whilst both had 
huge outbursts hundreds of years ago, P Cyg has been roughly 
constant in luminosity for many years, while r/ Car's variabil- 
ity is often explained by the expanding nebula moving across 



Fig. 10. Strength of depolarization against apparent variability over 
the last ten years. The symbols show the same as Figure 9. 



our line of sight, coupled with the effects of being in a highly 
eccentric binary system (Israelian & de Groot 1999; Davidson 
& Humphreys 1997). This leaves AG Car, HR Car and HD 
160529 as the only 'normal' Galactic LBVs. We may expect 
the MC LBVs to be a more homogeneous group, but as the data 
quality for the MC LBVs is not as good as for the Galactic ones, 
and because of the low sample numbers involved, we draw no 
conclusions as to differences between these two populations at 
the current stage. 

4.1.3. Wind asphericity and strong-amplitude 
variability 

A correlation can be found by looking at recent light-curves 
of the objects. Observations from the last - 15 years obtained 
from the AAVSO and the Tycho mission show that, of the ob- 
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jects showing definite line-effects, all show variability of > 1 
mag (see Fig. 10). Of the others, only S Dor shows this kind of 
variability, and this object was a borderline detection. We there- 
fore propose a link between wind asphericity and recent, strong 
variability. Figure 1 1 shows the light-curve of AG Car, going 
back 15 years. Marked on the plot are the H-D limit and the bi- 
stability jump (BSJ, see Sect. 1 and Lamers et al. 1995; Vink 
et al. 1999), under the assumption that the star varies at constant 
bolometric luminosity. It can be seen that AG Car has crossed 
its H-D limit and spent long periods very close to the BSJ dur- 
ing this time. The same is true for R 127 (Waagen 2005) and S 
Dor, although the assumption of constant bolometric luminos- 
ity may break down for S Dor (VG01). Whilst the luminosities 
of HR Car and R 1 10 mean they are below the H-D limit (leav- 
ing the uncertainties in their luminosities aside), both stars have 
shown recent strong variability (AAVSO, Waagen 2005), and 
are close to a second BSJ proposed by Lamers et al. (1995) at 
~ 10,000K, as Fe m recombines to Fe n. Proximity to the BSJs 
coupled with a variable effective temperature may play a part 
in producing wind asphericity. 

4.2. Interpretation of polarimetric variability 

We have shown that at least half of the LBVs observed dis- 
play evidence for intrinsic polarization. Of the 4 objects for 
which multi-epoch observations exist, all are polarimetrically 
variable, with both the degree of polarization and PA changing 
on timescales of ~ a week to months. As the bulk of the po- 
larization occurs within a couple of stellar radii, variable wind 
properties of an outflow moving at a few hundred km s _1 could 
produce variability on these timescales. Here we will discuss 
four different possible scenarios for producing such behaviour. 

4.2.1. Flip-flopping wind 

In this scenario, the wind flips from an equatorially-enhanced 
flow to a bi-polar flow or vice versa. Such a scenario could 
occur if the star was close to a BSJ and was rotating close 
to break-up, producing a temperature differential between the 
equator and the poles. As the star's global T e g changes, one of 
these regions may reach the BSJ temperature before the other, 
changing the geometry of the wind. As the plane of the equa- 
tor and poles are perpendicular to each other, at this point we 
may observe a 90° flip in intrinsic polarization. At intermediate 
stages we would expect the polarization of the two components 
to cancel each other out, and so we would observe lower levels 
of polarization but at PAs more or less aligned with one of the 
two axes. 

The results of SL94 showed that, for three separate obser- 
vations of AG Car over ~ 2 years, the PA of the intrinsic po- 
larization flipped roughly by 90° from 55° to 133° and back 
again. From Fig. 1 1 we can see that the star brightened by ~ 
0.2mags between their first and second observation, then had 
dimmed by ~ O.lmags by the time of their third. We can also 
see that the star was close to its BSJ at this time, assuming the 
star varies at constant luminosity. This alone is consistent with 
a flip-flopping wind. However, when one looks at the PAs mea- 



sured at 1994 (26°) and 2003 (84°), as well as those measured 
by Leitherer et al. (1994), these angles are completely unre- 
lated; whilst the degree of polarization remains in the range 
0.4-0.6%. This suggests that the results of SL94 may have been 
a result of limited temporal sampling, and that the flip-flopping 
wind cannot be the sole explanation of the polarization be- 
haviour. A similar situation is observed for HR Car and P Cyg, 
where the PA appears random with time and bears no relation 
to the level of polarization. 

4.2.2. Opacity changes within the wind 

As explained in Sect. 1, single scattering off a flattened wind 
would produce polarization perpendicular to the plane of scat- 
tering. However, if the opacity was to increase, multiple scat- 
tering effects would come into play. It has been shown by e.g. 
Angel (1969) and Wood et al. (1996) that multiple scattering 
in an optically thick, flattened wind can produce polarization 
parallel to the plane of scattering. If the density of the inner 
regions of the wind were to increase due to e.g. an increase 
in mass-loss rate, we may expect to see the PA of polarization 
change from perpendicular to parallel to the flattened wind's 
major axis. As with the flip-flopping wind, we would expect 
intermediate stages to have low degrees of polarization in one 
of the two planes due to cancellation effects. 

In Figure 1 1, the mass-loss rate for AG Car as determined 
by Stahl et al. (2001) are overplotted on the star's light curve. 
We can see that between the first two measurements by SL94, 
the mass-loss rate drops, coinciding with a 90° flip in the star's 
intrinsic polarization. However, at the time of SL94's third ob- 
servation, at which point the polarization had flipped back, the 
mass-loss rate had remained roughly constant. The seemingly 
unrelated PAs at similar levels of polarization of the later ob- 
servations are further reason to disregard this scenario for AG 
Car. The random temporal behaviour of HR Car and P Cyg also 
cannot be explained this way. 

4.2.3. Clumpy wind 

Here, the wind is not a smooth homogeneous outflow but in- 
stead consists of localized density enhancements, or 'clumps'. 
The bulk of the scattering would be due to the closest, dens- 
est clumps. If the clumps were ejected in a spherical distri- 
bution, or in a flattened wind which was observed at a low 
inclination, the PA of the intrinsic polarization would be re- 
lated to the projected angle between the star and the clump, 
and appear effectively random with time. Only in a flattened 
wind observed at a high inclination angle would the PA re- 
main roughly constant. Multiple clumps within the wind could 
serve to increase the intrinsic polarization or cancel each other 
out, making the degree of polarization also seemingly ran- 
dom with time. Wind-clumping has previously been invoked to 
explain the random polarimetric variations observed in Wolf- 
Rayet (WR) stars (Robert et al. 1989). 

AG Car, HR Car, and P Cyg are all described in this paper 
to have intrinsic polarization which appears to be essentially 
random with time, both in level and PA. This is consistent with 
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Fig. ll.The light curve of AG Car for 1989 
- 2004. Circles show data from the Tycho 
mission, crosses show validated data from 
the AAVSO records. The data has been re- 
binned to 50 days, and the error bars show 
the lcr standard deviations. The dates of 
our observation, the archive observation, 
plus the observations of SL94 have been 
marked. Also marked are the brightnesses 
at which the star is crossing the empirical 
H-D limit, and the bi-stability jump (BSJ) 
at ~ 21,000K (under the assumption that the 
star varies at constant bolometric luminos- 
ity). Also marked on the plot are the mass- 
loss rates derived by Stahl et al. (2001) (tri- 
angles). 



the clumpy wind scenario, with the clumps ejected in a roughly 
spherical distribution on the plane of the sky. Only R 127 has 
polarization with PA which remains roughly constant at 25 + 
10°. This is consistent with a clumpy wind ejected on a signifi- 
cantly flattened plane on the sky. The P Cygni profiles observed 
in the star's spectrum imply that material is being ejected to- 
wards us along the line-of-sight, therefore the favoured expla- 
nation of Schulte-Ladbeck et al. (1993) was that we are looking 
at an edge-on, equatorially-enhanced, clumpy wind. Our latest 
measurement supports this model. 



The strength of the polarimetric variability of LBVs is 
greater than that observed in WRs (Robert et al. 1989). It is 
interesting to note that Robert et al. find that the amplitude 
of the variability in WRs increases with decreasing terminal 
wind speed. This is consistent with the greater variability of 
LBVs (~ 1%, v ro ~ 200 km s" 1 ) compared with WRs (< 0.2%, 
Voo > 1000km s- 1 ). 



Investigations into the physical parameters of a single 
clump required to produce the observed levels of polariza- 
tion can be found in Nordsieck et al. (2001) and Rodrigues 
& Magalhaes (2000). The solid angle subtended by the clump 
must be large enough to scatter a significant fraction of the 
starlight. Nordsieck et al. (2001) find that the clump must 
be within ~ 2.5/? star and be ~ 20 times as dense as the am- 
bient wind; while Rodrigues & Magalhaes (2000) show that 
the clump must be a significant fraction of the size of the 
star (/? c iump > 0.5/? st ar)- Modelling the observed polarisation 
variablility of an LBV wind with multiple clumps at various 
distances in the wind and under a range of geometric distri- 
butions (inclination and opening angle) would yield valuable 
information about the overall wind geometry, as well as the 
clump production rate. 



4.2.4. Binarity 

If the variability of the intrinsic polarization were to be peri- 
odic, this would suggest that the polarimetric variations were 
due to interactions between the LBV and a binary companion, 
a la close WR-O systems. Here, the polarization is produced by 
the electrons in the WR star's fast, ionised wind scattering the 
light from the bright companion O star in the region between 
the two. The polarization describes a double loop in Q - U 
space about a central locus for each orbital cycle (see St. -Louis 
et al. 1987). The stars must be sufficiently close together in or- 
der for the WR's wind to intercept enough of the O star's light, 
and consequently the systems have periods of a few days to 
months. The temporal resolution of present polarimetric obser- 
vations is insufficient to detect such oscillations in the LBVs' 
polarization, and so this scenario warrants further invesigation. 
If a similar mechanism is responsible for the LBVs polariza- 
tion, the LBV would either have to play the role of the scatterer 
(i.e. be in a system with a hot star) or the illuminator (i.e. have 
a WR companion). 

If the LBVs were in a system with a WR star, we may ex- 
pect to see i) X-rays from the collision between the fast WR 
wind and the LBV wind, cf. the WC8 - evolved O star sys- 
tem HD 68273 (Schild et al. 2004); and ii) other periodic phe- 
nomena such as spectroscopic oscillations or the cyclical pho- 
tometric variability due to the comparable mass of the compan- 
ion. We note that of the four polarimetrically variable LBVs, 
only P Cyg has been detected as an X-ray source (Berghofer & 
Wendker 2000), and this has been rated as questionable by the 
authors. These stars have been subject to time-resolved pho- 
tometric and spectroscopic monitoring, but while oscillations 
have been detected, these are on characteristic timescales rather 
than strictly periodic, and are explained as stellar pulsations 
(see e.g. Stahl et al. 2001; Lamers et al. 1998). 

If the LBVs wind was electron-scattering the light of a 
companion, the luminosity of the companion would have to be 
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comparable to that of the LBV in order for the scattered light to 
make up the ~ 1 % polarization of the overall light we observe 
(the scattering of an O star's light in WR-0 binaries results in 
typically AP ~ 0.4%). It is unlikely that in this case the com- 
panion would not have been already detected. 

As an aside, it has been shown that after the periodic com- 
ponents have been subtracted, the polarimetric variations in 
WR-0 binary systems are no greater than for single WR stars 
(Robert et al. 1989). As these polarimetric variations are asso- 
ciated with wind-clumping in WR stars, this result implies that 
binarity does not play a significant role in the clumping of hot 
star winds. In short, binarity is an unlikely explanation for the 
observed polarimetric variability, either as a scattering mecha- 
nism or as a catalyst for wind-clumping. 

5. Conclusions 

From our observations, it is clear that LBVs show a high rate 
of wind-asphericity. We find evidence for intrinsic polarization 
in > 50% of those stars studied, compared with 25% and 20% 
for similar studies of O supergiants and Wolf-Rayet stars re- 
spectively. However, multiple observations reveal levels and 
angles of intrinsic polarization that are seemingly random with 
time. In investigating four possible explanations for this we ar- 
gue the evidence points away from the simple explanation of 
equatorially-enhanced or bi-polar flows. Instead, we interpret 
the variable intrinsic polarization as evidence for significant 
clumping within the wind. The random levels of polarization 
measured in R 127 are confined to position angles (PAs) of 
+ 10°, which can be explained as an equatorially enhanced, 
clumpy wind viewed edge-on. If viewed at lower inclination 
angles, the measured PAs may seem as random as those ob- 
served in AG Car, HR Car and P Cyg at limited temporal sam- 
pling. Therefore, time-resolved spectropolarimetric monitoring 
of LBVs would be required in order to detect axi-symmetry in 
their winds. 

We find that wind asphericity is more likely to be found in 
stars with strong Ha emission and strong recent variability (> 
lmag). As both these properties are linked to the stars' mass- 
loss rates they may also be linked to the clumping of the wind. 
With higher mass-loss rates, slower winds and lower effective 
gravities than their evolutionary neighbours, as well as residing 
in the bi-stability temperature regime, small changes in radius 
and/or effective temperature may lead to significant wind insta- 
bilities and inhomogeneities. This is consistent with the lower 
incidence of intrinsic polarization found in O supergiants and 
Wolf-Rayet stars. 
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Fig.A.l. A spectral atlas of early-type supergiants in the range 6240 - 6700 A. The prominent photospheric absorption lines of Sill and 
Hel, which appear to be good indicators of effective stellar temperature, have been identified. Spectral types of the stars were obtained from 
SIMBAD. The data in the atlas is available electronically through http : //vizier .u-strasbg . fr/viz-bin/VizieR. 
*HD 183143 is a candidate LBV. 



